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ABSTRACT 

The increasingly deep limit on the neutrino emission from gamma-ray bursts (GRBs) with IceCube 
observations has reached the level that could put useful constraints on the fireball properties. We first 
present a revised analytic calculation of the neutrino flux, which predicts a flux an order of magnitude 
lower than that obtained by the IceCube collaboration. For benchmark model parameters (e.g. the 
bulk Lorentz factor is T = 10 2 ' 5 , the observed variability time for long GRBs is t° h = 0.01s and the 
ratio between the energy in accelerated protons and in radiation is r/ p = 10 for every burst) in the 
standard internal shock scenario, the predicted neutrino flux from 215 bursts during the period of the 
40-string and 59-string configurations is found to be a factor of ~ 3 below the IceCube sensitivity. 
However, if we accept the recently found inherent relation between the bulk Lorentz factor and burst 
energy, the expected neutrino flux increases significantly and the spectral peak shifts to lower energy. 
In this case, the non-detection then implies that the baryon loading ratio should be r\ v < 10 if the 
variability time of long GRBs is fixed to t° h = 0.01s. Instead, if we relax the standard internal 
shock scenario but keep to assume r/ p = 10, the non-detection constrains the dissipation radius to 
be R > 4 x 10 12 cm assuming the same dissipation radius for every burst and benchmark parameters 
for fireballs. We also calculate the diffuse neutrino flux from GRBs for different luminosity functions 
existing in the literature. The expected flux exceeds the current IceCube limit for some luminosity 
functions, and thus the non-detection constrains r\ v < 10 in such cases when the variability time of 
long GRBs is fixed to t° h = 0.01s. 
Subject headings: gamma-ray bursts — neutrino 



1. INTRODUCTION 

GRBs have been proposed as one of the potential 
sources for ultra-high energy cosmic rays (UHECRs) with 
energy up to > 10 20 eV (e.g. Waxman 1995; Vietri 
1995; Waxman & Bahcall 2000; Dai & Lu 2001; Der- 
nier 2002; Murase et al. 2006), given the hypothesis 
that the fireball composition is proton-dominated and 
the protons get accelerated in the dissipative fireballs. 
Interactions of protons with fireball photons will pro- 
duce a burst of neutrinos with energies of ^PeV (e.g. 
Waxman &Bahcall 1997; Guetta et al. 2004; Dermer 
& Atoyan 2006), so detection of such neutrinos would 
prove the presence of cosmic ray protons in the fireball. 
Despite that large progresses in the studies of GRBs 
and their afterglow have been made recently, the com- 
position of the jet, whether it is proton-electron domi- 
nated or Poynting-flux dominated, is largely unknown. 
Both baryon-dominated fireball shock model (e.g. Rees 
& Meszaros 1994; Paczyhski & Xu 1994) and magnetic 
dissipation model (e.g. Narayan & Kumar 2009; Zhang 
& Yan 2010) have been proposed for the central engine 
of GRBs. In the magnetically dominated outflow model 
for GRBs, the non-thermal proton energy fraction may 
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be low (but see Giannios 2010), while in the baryon- 
dominated outflow model, it is natural to expect proton 
acceleration via shock dissipation of the kinetic energy. 
Since the flux of neutrinos depends on the energy frac- 
tion of protons in the fireball for a given burst energy, the 
flux or limit of the neutrino emission could constrain the 
proton energy fraction and in principal provide a useful 
probe of the jet composition. The proton energy frac- 
tion is also crucial to know whether GRBs could provide 
sufficient flux for UHECRs. 

The Kilometer-scale IceCube detector is the most sen- 
sitive neutrino telescope in operation, although no pos- 
itive neutrino signal has been detected so far (Abbasi 
et al. 2010, 2011a; The IceCube collaboration 2011). 
The IceCube operations with the 22, 40 and 59 strings 
configurations all yield negative results, which have put 
more and more stringent constraints on neutrino emis- 
sion from GRBs. The analysis were performed for both 
point source search from individual GRBs and diffuse 
emission from aggregated GRBs (Abbasi et al. 2011b). 
According to the IceCube collaboration, with 40-string 
configuration operation between 2008 and 2009, IceCube 
reaches a sensitivity at the level of the expect flux from 
GRBs and the combined upper limit of IceCube 40-string 
and IceCube 59-string analysis is 0.22 times the expected 
flux (Abbasi et al. 2011a). Based on this, the IceCube 
collaboration argued that the UHECR-GRB connection 
is challenged (Abbasi et al. 2011a). 

The calculation by the IceCube collaboration (ICC 
hereafter) is based on the formula in the appendix of the 
paper Abbasi et al. (2010) and benchmark parameters 
for the internal shock model of GRBs. However, as also 
pointed out by Li (2011) and Hummer et al. (2011), we 
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will show that the normalization procedure used by ICC 
overestimates the neutrino flux. In calculating the pho- 
ton number density, ICC also approximates the energy 
of all photons by the break energy of the photon spec- 
trum, which originates from the assumption in Guctta 
et al. 2004 (GT2004 hereafter). To get a more accurate 
estimate of the expected neutrino flux, we first present 
a refined analytic calculation in § 2.1, revising the above 
approximations, and perform a numerical calculation in 
§ 2.2 taking into account three main energy loss chan- 
nels for protons interacting with burst photons. In § 3, 
we study the effects of non benchmark parameters on the 
neutrino emission, such as the dissipation radius and the 
bulk Lorentz factor of the fireball. In § 4, we calculate the 
accumulative diffuse neutrino emission from GRBs and 
confront it with the IceCube limit on diffuse neutrinos. 
At the end, we give our conclusions and discussions. 

2. GRB NEUTRINO SPECTRA 

Based on the assumption that protons and electrons 
are accelerated in the same region of a GRB, pro- 
tons interact with photons emitted by electron syn- 
chrotron emission or inverse- Compton emission predom- 
inantly produce the charged and neutral pions. The 
charged pion subsequently decays to produce 4 final 
state leptons, via the processes ir^ — > u^D^)^ — > 
zy /i(f'/i)e + (e _ )^ e (i/ e )P Al (f AI ), which approximately share 
the pion energy equally. Denoting 5ft as the ratio be- 
tween the charged pion number to the total pion number, 
the fraction of the proton energy lost into each lepton is 
§/ P7 , where / P7 is the fraction of the protons energy lost 
into pions. For the proton with energy of e p — 7 p 77i p c 2 , 
the photomeson interaction timescale can be calculated 
by (Waxman & Bahcall 1997) 
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where er p7 (e 7 ) is the cross section of photomeson inter- 
action for a photon with energy e 7 in the proton-rest 
frame, £(e 7 ) is the inelasticity, e 7) th is the threshold of 
photon energy and -p 1 is the GRB photon spectrum in 

fluid- rest frame (Waxman & Bahcall 1997). The fraction 
of protons energy loss into pions is 



fp'y 1 CXp( tdyn/tpj*) : 
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where idyn = R/(Tc) is the dynamic timescale. Just for 
simplicity, we use idyn as the interaction time. Gener- 
ally speaking, proton cooling timescales can be shorter 
than the dynamical timescale (see Murase & Nagataki 
2006a, b, for details), but we do not consider such com- 
plicated effects for the purpose of testing the standard 
model suggested by Waxman & Bahcall (1997) (Hum- 
mer et al. 2012). 

Approximating that neutrinos produced via photome- 
son interaction by proton with energy e p have constant 
energy, the spectrum of neutrinos from the decay of sec- 
ondary particles can be obtained by, without considering 
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the oscillation, 
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where -4^ is the spectrum of protons and 5ft is the ra- 
tio between the charged pion number to the total pion 
number as defined above equation {l}. The subscript i 
represents different flavors of neutrinos, i.e., i — f M , v^, v e 
for muon neutrinos produced via the decay of pions, 
antiunion neutrinos P„ and electron neutrinos v e pro- 
duced via the decay of muons, respectively. If the cool- 
ing timescale of pions or muons is smaller than their 
lifetime, pions or muons have the probability to cool 
down before decay. Then the neutrino flux will be sup- 
pressed by a factor of £ ff = 1 — cxp (— ^, syn /r w ), where 
^, sy n = 3.3 x 10~ 3 sL~ g 2 r| 5 i?f 4 e~ is the synchrotron 

cooling timescale, and tv = 2.6 x 10 _8 S7, r = 186se,r ; EeV 
is the lifetime of pions, whose energy is e„ = 0.2e p . Here 
we assume the fraction of electron energy and magnetic 
field energy are same, i.e., e e — es- Similarly, the sup- 
pression due to muon cooling is = 1 — exp (— t tl , syn /T tl ) , 

where t^ syn = 1.1 x 10 _3 s J L~5 2 r| 5 i?? 4 e^ 1 EoV and r M = 

2.1 x 10 4 se M ,EcV with the muon energy e M = 0.15e p . 
The suppression factor 9e{e p ) is a combination of Cr( e p) 
and Cm( £ p)> i- e -j ^( e p) = Cr(e P ) and ^(„ e )(e p ) = 

Considering the neutrino oscillation effect, the spec- 
trum of muon neutrinos (including antimuon neutrinos) 
detected on Earth is approximated as ( Nagataki et al. 
2003; Particle Data Group 2004; Kashti & Waxman 2005; 
Murase 2007; Li 2011) 
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where are produced via the decay of secondary pions, 
and v e and are produced via the decay of secondary 
muons. Therefore, by inserting equation ([3]) into equa- 
tion Q, one can calculate the total spectrum of muon 
(including antimuon) neutrinos detected on Earth nu- 
merically via the following equation: 



dn u Jftjep) dn p 

e v -j^-ae u — — - — J pl {e p )y u {e p )e p -^-ae p , 
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with the factor 9 u {e p ) = 0A(^ (e p ) + 0.6Q(e p )^ t _ l (e p ) ac- 
counting for the neutrino oscillation and the cooling of 
secondary particles. 
The GRB photon distribution can be described by 
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where e 7; fc is the break energy (in the fluid-rest frame) 
of the photon spectrum, q — a for e 7 < e 7 j, and 
q = (3 for e 7 > e 7 b- The normalized coefficient is 
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for p ^ 2, where e 7jm i n and e 7jmax are the minimum 
and maximum energy of the photons. For /3 = 2, 

yi = (l (^y^ 2 ) /(-o + 2) + ln(^ f ). Here- 

after, we adopt assumptions that e 7 m j n = IkeV and 
e 7jmax = lOMeV as ICC did (Abassi et al. 2010). If we 
know redshift of bursts, the correction should be properly 
taken into account in estimating the luminosity. 

The proton number per energy interval can be de- 
scribed by = N p e~ s with s being the power-law index 

and N p being the normalized coefficient. The normalized 
coefficient of the injected proton spectrum can be calcu- 
lated by N p = Epj y*'2Z dt P e l~ S = E vl ln T~ ( we set 
s = 2 hereafter, as predicted by the shock acceleration 
theory) with E p being the total energy in protons, e Pj: 



and 



being the minimum and maximum energy of 



accelerated protons, respectively. We introduce a factor 
rj p , denoting the ratio of the energy in accelerated pro- 
tons to the radiation energy, then the total proton energy 
E p is (Murase & Nagataki 2006a) 



E p — r] p E ie 
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where .Ei SO is the isotropic energy of the burst, which 
is obtained from the observed fluence F° h in the en- 
ergy band of 1 keV - 10 MeV and the redshift of the 
burst. According to Waxman (1995) and Li (2011), 
for mildly-relativistic GRB internal shocks, we assume 
the minimum proton energy as £p b min — TmpC 2 = 
3.0 x 10 n r2.5/(l + z)eV , and the maximum proton 
energy due to synchrotron cooling is £p b max = 4.0 x 

ia20t^3/2 r.l/2 1/4 -1/4 -1/2 T -1/4 , / 1 . \ \ 7 ^1 

10 r 2 :5#i4 £ e e B 9i L 7 ,52 /(l+«)eV, with gi > 1 
being a factor accounting for the uncertainty in the par- 
ticle acceleration time. 



2.1. Analytical Calculation 
2.1.1. Neutrino spectrum in the general dissipation scenario 

In this subsection, we treat the dissipation radius as 
a free parameter, since the exact dissipation mechanism 
of GRBs is not established. There are suggestions that, 
besides the standard internal shock model, the prompt 
emission arises from the dissipative photosphere or arises 
at much larger radii where the magnetic-dominated out- 
flow is dissipated through magnetic dissipation processes, 
such as reconnection (e.g. Narayan & Kumar 2009, Ku- 
mar & Narayan 2009, Zhang & Yan 2010). 

For the analytical calculation, we adopt the A reso- 
nance approximation as in Waxman & Bahcall (1997) 
and Guetta et al. (2004), where the cross section peaks 
at the photon energy e 7 ~ e pC ak = 0.3GeV in the proton- 
rest frame. If i^yn < Vy ^ ne conversion fraction is 
approximated as / P7 ~ fdyn/Vy = R/fictp^). Adopt- 
ing the A resonance approximation (Waxman & Bahcall 
1997), the fraction of proton energy converted into pion 
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: ATrripC 2 or e° b min ~ WFmpC 2 



are adopted, since the relative Lorentz factor is order of 1-10, where 
higher values favor efficiency internal shocks. 
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Note 



that the above approximation is valid when 
the radius is not too small, i.e., R > 1.1 x 

10 13 (£) (sk) (i^) W^V^Cm. 

As t v — 0.0 5e p for the A resonance approximation, 
from equation (|10p . we can obtain the break energy of 
neutrino spectrum corresponding to the photon spectral 
break 



e°% = 7.5 x 10 5 GeV(l + zy 2 T 2 2 5 e 
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The cutoff energy of muon neutrino spectrum due to the 
pion cooling can be obtained by setting i w , sy n = TV, 



ob = 3 : 3xl0 8 r _V2 T . 2 
(1 + 



V52 risfluGeV. 
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For antimuon neutrinos (and electron neutrinos) pro- 
duced via the decay of muons, an extra break is caused 
by the muon cooling, which is 



ob 
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where the subscript A represents antimuon neutrino P M 
or electron neutrino v e . 

Assuming that the fraction of the amount of charged 
pions is 5ft = 1/2, from equations ([3|) and ((9]), we get 
the spectrum of muon neutrinos produced by the pion 
decay, 
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Hereafter, for brevity, we abandoned the coefficient fei and k t 
in the following equations, since k± ~ 1 and kg ~ approximately, 
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where k( 
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with £ represents three flavors of neu- 



trinos, i.e., electron neutrinos and muon(antimuon) neutrinos. But 
we still adopt that in our analytic calculations. 
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Similarly, the spectrum of antimuon (and electron) neu- 
trinos produced by muon decay is approximated by 
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.Then, one can obtain the + spectrum after consid- 
ering the neutrino oscillation effect by substituting equa- 
tions (HU) and CE5]) into equation (gj). 



and the spectrum of antimuon(electron) neutrinos pro- 
duced via muon decay is approximated by 
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where the cutoff energies are 
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with A representing antimuon and electron neutrinos [y^ 
and v e ) produced by muon decay, and 



2.1.2. Neutrino spectrum in the internal shock scenario 

In this subsection, we assume the standard inter- 
nal shock scenario with the dissipation radius at R = 
2T 2 ct° b /(l -I- z), where t° h is the observed variability 
timescale of GRB emission. The conversion fraction / P7 
is given by 
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Then the spectrum of muon neutrinos produced via pion 
decay is approximated by 
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Substituting equations (fTT]) and (fl8|) into equation d4} , 
we can obtain the neutrino spectrum analytically. To il- 
lustrate the difference between our calculation and the 
ICC calculation, we calculate the neutrino spectrum for 
one typical GRB with benchmark parameters, shown in 
Figure 1. Compared with the ICC calculation (the dark 
gray solid line) , our spectrum (the purple solid line) con- 
sists of more structures resulting from the sum of the con- 
tributions by three flavors of neutrinos, for which both 
pion cooling, muon cooling and the oscillation effect are 
considered. Furthermore, the flux level predicted by our 
modified analytical calculation is a factor of ~ 20 lower 
than that obtained by ICC (Abbasi et al., 2011). This 
mainly arises from two differences in the calculation: 

(I) We use equation ©, where the conversion fraction 
/ P7 is a function of the proton energy e p as shown by 
equation (I16p . to normalize the neutrino flux to the pro- 
ton flux, which means that only a fraction of protons can 
produce neutrinos efficiently. This corrects ICC's inac- 
curate use of energy-independent conversion fraction in 
the normalization of the neutrino flux (Li 2011; Hummer 
et al. 2011; Murase et al. 2012). The calculation of 
Guetta et al. (2004) Q normalized the flux based on the 

7 Guetta et al. (2004) calculated the neutrino spectrum by as- 
suming a flat high-energy electron spectrum (i.e. dN e /d'y £ oc 7<T ) 
and using an electron equipartition fraction e£F that represents the 
ratio of the nonthermal electron energy over one energy decade to 
the UHECR energy over one energy decade, and the neutrino flux 
is normalized by eldN„/de v = (l/8)(l/e^)(F° b /lnlO)/^ (see their 
Eq. A19). Note that other normalization procedures are also pos- 
sible, and this e£? is typically larger than the conventional e e that 
is defined as the ratio of the total nonthermal electron energy to 
the total internal energy (including both thermal and nonthermal 
protons) . 
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differential spectrum so that it does not suffer from this 
problem. The spectrum obtained with the calculation of 
Guetta et al. (20040 is shown by the blue solid line in 
Figure 1, where the middle of equation (A. 15) in their 
paper is used assuming the bolometric luminosity as the 
luminosity at the break energy. The flux is lower than 
the ICC result (the dark gray solid line in Figure 1) by 
a factor of ~ 4. 

(II) In calculating the photon number density, we con- 
sider the photon energy distribution according to the real 
photon spectrum, reflected by the normalized coefficient 
Ay ~ U 1 , where y\ is shown in equation ([7]), while 

ICC approximate the energy of all photons by the break 
energy of the photon spectrum. This leads to a flux a fac- 
tor of ^3-6 lower than the ICC result for typical a and 
f3 values and e 7 b ~ 100 — lOOOkeV in the GRB spectrum. 

Finally, we get a neutrino spectral flux (purple solid 
line) lower than the one predicted by ICC (dark gray 
solid line) by a factor of ~ 20 for a typical GRB, as 
shown in Figure 1. Of course, the suppression factor is 
different for GRBs with different parameters. 

2.2. Numerical Results 

Besides the baryon resonance, the direct pion process, 
multi-pion process and the diffractive scattering also con- 
tribute to the total P7 cross section. Therefore, in our nu- 
merical calculation, we adopt a more precise cross section 
for photomeson interaction, including three main chan- 
nels, i.e., the A resonance process, direct pion process 
and multi-pion process. For simplicity, we assume that 
the inelasticity is £ = 0.2 for e 7 < 983MeV and £ = 0.6 
for e 7 > 983MeV(Atoyan & Dermer 2001). Then, insert- 
ing the cross sections and inelasticity into equation (fTJ 
and using equation @, we can get the fractions of pro- 
tons energy converted to pions. The average fractions 
of the charged pions are set to 3?a = 1/3 for A Res- 
onance process, SRdir = 2/3 for direct pion production, 
and 5R mu i = 2/3 for multipion production, based on nu- 
merical investigations (Mucke et al. 1999, 2000, Murasc 
& Nagataki 2006a, b, Murase et al., 2006, Murase, 2007, 
Baerwald et al. 2011). In addition, we assume the neu- 
trino energy is e„ = 0.05e p for A-resonance and direct 
pion production channels, and e„ = 0.03e p for multipion 
production channel. Inserting these quantities into equa- 
tions ((J), @ and then ([5]), we can obtain the spectrum 
of neutrino emission produced via the three dominant 
channels. Although this simplified numerical approach 

8 In our paper, we used the nonthermal baryon-loading pa- 
rameter rip, defined as the ratio between the total energy in 
accelerating protons and the total radiation energy at observed 
bands, to normalize the neutrino flux through e^,dN v /de v = 
(l/8)(j?p/ln(e p , max /e Pimin ))_F° f, /p 7 . The lines of "modified Guetta 
et al. (2004) " in Figures 1 and 2 are obtained by using 
the benchmark value of r\ v = 10, which corresponds to = 
ln(ep, m ax/e p ,min)/hi(10)(l/?7p) ~ 1 in Eq.(A19) in Guetta et 
al. (2004). Choosing = 0.1 in Guetta et al. (2004) 
would correspond to a higher baryon loading factor of r\ v = 

ln(ep jm ax/ep,min)/ln(10)(l/e5?) ~ 100 in our paper. Note that the 
GRB-UHECR hypothesis might suggest such high loading factors 
if the local GRB rate is 0.01 — O.lGpc - 3 yr — 1 or the local lumi- 
nosity density is significantly smaller than ~ 10 44 ergMpc _3 yr — x , 
and those relatively optimistic cases are constrained by IceCube 
under reasonable assumptions, as expected in Murase & Nagataki 
(2006b)." 



is different from more detailed, fully numerical calcula- 
tions (Murase & Nagataki 2006a,b; Murase 2007; Baer- 
ward et al. 2011), it is enough for our purpose and saves 
calculation time. 

In Figure 1, we also show this numerical result (the red 
solid line) for comparison. The flux obtained with this 
numerical calculation is about 2-3 times larger than the 
analytical result in the energy range 10 5 GeV-3 x 10 6 GeV 
for typical GRB , consistent with Murase & Nagataki 
(2006a, b) and Baerwald et al. (2011). So the analytic 
calculation presented in the above section can still be 
used as a rough approximation. 

2.3. Confronting the calculations with IC59+40 
observations 

215 GRBs are observed during the operations in the 40- 
string and 59-string configurations of IceCube, yielding 
negative results. In this section, we calculate the neu- 
trino flux for the same 215 GRBs using the same burst 
parameters as ICC. The information for these samples 
was taken from the website grbweb. icecube. wise. edr0 and 
GCN. For the unmeasured parameters, we adopt the 
same average values as in Abbasi et al. (2010, 2011a). 
We assume a ratio of proton energy to radiation en- 
ergy r/p = 10, equivalent to l/f e = 10 in ICC calcu- 
lation. We adopt the internal shock model with shock 
radius at R = 2T 2 ct° h /(l + z), where the Lorentz fac- 
tor is r = 10 2 ' 5 and the observed variability timescale 
is t° h = 0.01s for each long GRB and t° h = 0.001s for 
each short GRB as the ICC did. The diffuse neutrino 
flux can be obtained by the total neutrino fluence for 
215 individual GRBs timing a factor 7.83 x 10 _9 sr _1 s _1 
with the assumption that a total of 667 uniformly GRBs 
are generated per year. By adopting the effective area 
of IC59 and IC40 as a function of zenith angle, we can 
calculate the expected number of neutrinos with energy 
from 10 5 GeV to 3 x 10 6 GeV (Abbasi et al. 2010, 2011a). 
The corresponding combined 90% confidence level (CL) 
upper limit spectrum can be obtained by assuming that 
the limit amount of neutrino events is nn m ~ 1.9 in the 
energy range from 10 5 GeV to 3 x 10 6 GeV as in the 
IceCube collaboration (2011). 

In Figure 2, we show the 215 neutrino spectra (light red 
thin solid lines) for individual GRBs. The sum of them 
is presented as the red thick solid line, which is about 
a factor of ~ 10 lower than that predicted by ICC (the 
dark gray solid line; The IceCube collaboration 2011). 
As we can see, the neutrinos predicted by ICC2010 and 
GT2004 (dark gray solid line and blue solid line) are 
above the 90% CL upper limits (dark gray dashed line 
and blue dashed line) for the combined IC40 and IC59 
data analysis, while our predicted neutrino flux is below 
the corresponding upper limit (the red dashed line). We 
find that the total expected number of neutrinos with 
energy from 10 5 GeV to 3 x 10 6 GeV is 0.74, which is 
36% of the 90% CL upper limit. 

3. CONSTRAINTS ON THE FIREBALL PROPERTIES 

3.1. Uncertainty in the dissipation radius 

In § 2.1.1, we do not assume a specific dissipation 
model for GRB emission, but leave the dissipation radius 

9 For some information from this website with typo, which is 
different from GCN, we take the information from GCN. 
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as a free parameter. Internal shocks happen over a wide 
parameter range (e.g. Nakar & Piran 2002), and even 
larger radii are suggested in some other dissipation mod- 
els (e.g., Narayan & Kumar 2009, Zhang & Yan 2010). 
In Figure 3, we show the total neutrino spectra for 215 
GRBs by assuming some fixed dissipation radii for every 
GRBs in the range R = 10 12 cm - 10 16 cmH- It shows 
that the neutrino flux for the case of R = 10 12 cm (the 
black solid line) would exceed the corresponding IceCube 
upper limit (the black dashed line) as long as the baryon 
loading factor is sufficiently greater than unity. If we 
fix rjp = 10, the non-detection requires that the dissipa- 
tion radius should be larger than 4 x 10 12 cm. We note 
that, when the emission radius is too small, the max- 
imum energy of accelerating particles is limited due to 
strong photohadronic and/or radiation cooling, and neu- 
trino emission can be more complicated due to strong 
pion/muon cooling, so a more careful study is needed 
to obtain quantitative constrains on rj p in this regime. 
On the other hand, the larger dissipation radius leads 
to lower neutrino flux and higher cooling break energy 
according to equations (12) and (13). The shift of the 
first break to higher energies for larger dissipation radii 
is due to that those GRBs with a > 1, whose neutrino 
spectral peaks locate at the cooling breaks, contribute 
dominantly to the neutrino flux. 

3.2. Uncertainty in the bulk Lorentz factor 

In the previous subsections, we took either the variabil- 
ity or dissipation radius as a principal parameter, given 
a Lorentz factor, i.e. V — 10 2 5 . For those bursts with- 
out measured redshift, we took L 1 = 10 52 ergs _1 for the 
peak luminosity, as done by ICC. However, it was found 
recently that, the bulk Lorentz factor could vary signif- 
icantly among bursts, and there is an inherent relation 
between the Lorentz factor and the isotropic energy or 
the peak luminosity (Liang et al. 2011; Ghirlanda et. al. 
2011). As shown by Equations CLZ]) and (JTHJ), the neu- 
trino flux is very sensitive to the bulk Lorentz factor, so 
we can use the inherent relation to obtain more realistic 
values for the Lorentz factors and, hence, more reliable 
estimate of the neutrino flux. 

By identifying the onset time of the forward shock from 
the optical afterglow observations, Liang et al. (2011) 
and Lv et al. (2011) obtain the bulk Lorentz factors 
for a sample of GRBs. They further found a correlation 
between the bulk Lorentz factor and the isotropic energy 
of the burst, given b}0 

T L = 118£° so 26 52 . (22) 

Ghirlanda et. al. (2011) revisit this problem with a large 
sample and obtain a relation as 

T G - 29.8£? so 51 52 . (23) 

Compared with the benchmark model which assumes 
F = 10 2 ' 5 for all bursts, the value of T obtained from 

10 We note that, if the radius is smaller than the photosphere 
radius, the neutrino emission produced by p — p interactions be- 
comes important (Wang & Dai 2009; Murase 2008), which is not 
considered here. 

1 We adopt only the center value for the relationships presented 
hereafter. 



these relations is lower for bursts with the isotropic en- 
ergy E iso < (4.4 - 9.4) x 10 53 erg. 

Ghirlanda et. al. (2011) also obtained the relation 
between the bulk Lorentz factor and peak luminosity, 

r Gi =72.1L°$. (24) 

The value of T obtained from this relation is lower than 
10 2 - 5 for bursts with luminosity L p < 2.0 x 10 53 ergs~ 1 . 

Yonetoku et al. (2004) and Ghirlanda et al. (2011) 
found an inherent relation between the peak energy of 
photon spectrum and the peak luminosity. Therefore one 
can obtain the peak luminosity from the observed break 
energy of the photon spectrum e° b b and the redshift z by 

adopting the e° b b — L 7 relation, which is 

i 7 G,52 = 7.54[ £ ° b biMcV (l + z)] 1 - 75 , (25) 

derived by Ghirlanda et. al. (2011). 

We use the above inherent relations to calculate both 
the Lorentz factor T and the peak luminosity i 7 , and 
then calculate the neutrino flux produced by the same 
215 GRBs, which is shown in Figure 4. The main dif- 
ferences in the neutrino spectrum resulted from using 
different choices of the value of the bulk Lorentz factor 
can be summarized as: 

(i) The peak energy of the neutrino spectrum shifts 
to lower energy for the models adopting the relations in 
Ghirlanda et al. (2011) and Lv et al. (2011). This is 
due to that for the majority in the 215 GRBs, the values 
of r derived with these inherent relations are lower than 
the benchmark value 10 2 5 , which leads to a lower peak 
energy in the neutrino spectrum according to equation 
(jTTJ) . Also, the cutoff energy shifts to lower energies for 
these models according to equations (fT"9| and (|20|) . 

(ii) The peak flux of the neutrino spectrum increases 
for the two models which adopt the inherent relations. 
This is due to that, a lower Lorentz factor leads to 
higher neutrino production efficiency in the internal 
shock model. The predicted neutrino flux for both mod- 
els adopting Ei so — T relations in Ghirlanda et al. (2011) 
and Lv et al. (2011) exceed the IceCube upper limit, 
which implies a baryon ratio r\ v < 10 if t° b = 0.01s for 
long GRBs is correct. 

In the left panel of Figure 4, we take the redshift 
z = 2.15 for those long GRBs without measured red- 
shift, the amount of which is about 84% of the total 
amount of GRBs. For the benchmark model, the neu- 
trino flux would not be affected significantly since it is 
independent of the redshift according to equations (jTTJ) 
and (|18p. However, the value of redshift can affect the 
flux of the neutrino spectrum for models adopting the 
inherent relations. For a fixed observed fluence of the 
7— ray emission, a smaller redshift will lead to a smaller 
peak luminosity or isotropic energy. As a result, the 
Lorentz factor, derived from the inherent relations, will 
be lower, which leads to a higher neutrino production ef- 
ficiency. As shown in the right panel of Figure 4, the flux 
for the two models adopting inherent relations increase 
if we take z = 1 for those long GRBs without measured 
redshifts. 

In the above discussions, we have implicitly assumed 
the baryon ratio r) p = 10 for GRBs, in accordance with 
the notation / e =0.1 in Abbasi et al. (2010). This value 
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comes from the assumption that the radiation efficiency 
for GRBs is typically 0.1 and that most of the dissipated 
internal energy goes into the accelerated protons. How- 
ever, the fraction of energy in protons, and hence the 
value of r]p, are not well-known. The null result of Ice- 
Cube observations allows us to put some constraints on 
this value. Define r) p ^ c as the critic value, above which 
GRBs would be detected by the corresponding IceCube 
configurations. In Table 1, we list the corresponding 
value of r/p^c for the combined IC40 and IC59 analysis. 
However, one should keep in mind that rj PtC depends on 
the choice of t° h , so it should be larger for larger values 

4. DIFFUSE NEUTRINO EMISSION FROM GRBS 

Recently, IceCube also reported observations of dif- 
fuse neutrinos by the 40-string configuration. The 
non-detection yields an upper limit of 8.9xlO _9 GeV 
cm s sr , for the diffuse neutrino flux assuming an 
E~ 2 neutrino spectrum (Abbasi et al. 2011b). The ex- 
pected diffuse GRB neutrino flux can be obtained by 
summing the contributions of all GRBs in the whole uni- 
verse. To this aim, we also take into account the number 
distribution of GRBs over the luminosity (i.e., the lu- 
minosity function) as well as the number distribution at 
different redshifts. Although it is not relevant if the GRB 
parameters (T and R) do not depend on the luminosity, 
it affects results when one adopts specific relations such 
as equations ([22]) or ([23]) or (|2B]) . 

We employ three different luminosity functions and the 
corresponding source density evolution functions to de- 
scribe the distribution of the GRB number over luminos- 
ity and redshift. One luminosity function is suggested by 
Liang et al. (2007, hereafter LF-L), 



dL 1 



L 7 b 



(26) 



where po = 1.2Gpc _3 yr _1 is the local event rate of 
GRBs, and $o is a normalization constant to assure the 
integral over the luminosity function being equal to the 
local event rate po- This luminosity function breaks at 
L 7 b = 2.25 x 10 52 ergs~ 1 , with indices ai = 0.65 and 
o?2 = 2.3 for each segment. The normalized number 
distribution of GRBs with redshift used by Liang et al. 
(2007) in obtaining this luminosity function is (Porciani 
&Madau 2001) 



S(z) = 23- 



„3.42 



,3.42 



22.0 



(27) 



Wanderman & Piran (2009) also suggested a luminos- 
ity function in the form of a broken power-law (hereafter 
LF-W) 



dN _ { (t^)-" 1 L < L 7b , 
5^- P0 * O i(^)- L>L lh , 



(28) 



where po = 1.3Gpc~ 3 yr _1 , ol\ = 1.2, a 2 = 2.4 and break 
luminosity L 7 b = 10 52 5 ergs _1 . The corresponding nor- 
malized number distribution with redshift is described 

by 

(l + z) 2 - 1 z<3, 
(1 + z)- 1 - 4 z>3. 



S(z) 



(29) 



Another luminosity function we consider here is sug- 
gested by Guetta & Piran (2007, hereafter LF-G), which 
is in the same form as that of Wanderman & Pi- 
ran (2009), but with different parameters, i.e., po — 
0.27Gpc" 3 yr"\ a x = -1.1, a 2 = -3.0 and L 7b = 
2.3 x 10 51 ergs _1 . It implies a smaller local event rate 
and fewer GRBs at the high luminosity end. This lumi- 
nosity function is obtained based on the assumption that 
the rate of GRBs follows the star formation history given 
by Rowan- Robinson (1999), i.e. 



S(z) = 



10 
10 



,0.752 
,0.75 



Z < 1, 
Z>1. 



(30) 



Denoting the differential neutrino number generated 
by a GRB with luminosity L 1 at local redshift z by 
dn v /de lll the injection rate of neutrinos per unit time 
per comoving volume then can be obtained by 



' l " (Ly^^^-iL^dLy, (31) 



where p(z) = p S(z) is the event rate density in the rest 
frame. Considering the cosmological time dilation and 
the particle number conservation, a neutrino with energy 
e v observed at the Earth must be produced at redshift z 
with energy (1 + z)e u and *„(e° b )de° b = (1 + z)^ v [(l + 
z)e° h ]de° h . The total observed diffuse neutrino flux can 
then be integrated over redshift, 



dN b 



I^ 1 



cdz' 



(32) 



where H(z) = H /y/(l + z) 3 fl M + is the Hubble 
constant at redshift z. Here we set i 7 min = 10 50 ergs~ 1 , 
£ 7 ,max = 10 54 ergs _1 and z max = 8. 

The results of diffuse neutrino flux are shown in Figure 
5. In this plot, we show the diffuse muon neutrino spec- 
tra for both the case in which the bulk Lorentz factor 
of the GRBs follows the inherent relation suggested by 
Ghirlanda et al. (2011) (solid lines), as described in Sec. 
3, and the case in which a constant value of Lorentz fac- 
tor (F = 10 2 5 ) is assumed for all GRBs (dashed lines). 
The photon spectrum of GRBs is assumed to be a bro- 
ken power-law spectrum described by equation ([5]) with 
the a = 1 and (3 = 2, and the break energy of photon 
spectrum is calculated from the peak luminosity of GRBs 
via the relationship shown in equation (|25[) . One can see 
from this plot that using different luminosity functions 
and associated source density evolution functions leads 
to very different flux of diffuse neutrinos. The expected 
flux for the LF-L function slightly exceeds the IC40 upper 
limit, while for the luminosity functions of LF-W and LF- 
G, the predicted flux is undetectable even with one-year 
full operation of IceCube. Particularly, the LF-G func- 
tion results in a very low flux. This is not only because 
that the local event rate for this luminosity function is 
much lower, but also because that much more GRBs lo- 
cate at the low luminosity end, which contribute lower 
neutrino flux than more luminous ones. To require the 
diffuse neutrino flux not to exceed the IC40 upper limit, 
we have rj p <8 for the luminosity function of LF-L. For 
results in cases where the GRB parameters do not de- 
pend on the luminosity, see Murase & Nagataki (2006a), 
Murase et al. (2006), and Gupta & Zhang (2007). 



5. CONCLUSIONS & DISCUSSIONS 

The non-detection by the increasingly sensitive detec- 
tor IceCube has provided interesting implications for var- 
ious theoretical predictions of neutrino emission from 
GRBs. The IceCube collaboration reported that the Ice- 
Cube 40-string and 59-string configurations have reached 
the sensitivity below the theoretical expectation, which, 
if true, would challenge the view that GRBs could be the 
sources for UHECRs. However, as also shown by previ- 
ous works, we show that the IceCube collaboration used 
an overestimated theoretical flux in comparison with the 
IceCube instrument limit. We therefore revisit the ana- 
lytic calculation of the neutrino flux, considering the re- 
alistic photon energy distribution in calculating the num- 
ber density of fireball photons (instead of using the bolo- 
metric luminosity as the luminosity at the break energy) , 
and using the appropriate normalization for the proton 
flux to evaluate the neutrino flux. 

Using the modified formulas, we calculate the expected 
neutrino flux from the 215 GRBs observed during the op- 
erations of IceCube 40 and 59 strings configurations, as- 
suming the same benchmark parameters as that used by 
the IceCube collaboration. The flux is about 36% of the 
90% CL upper limit, consistent with the non-detection 
of IceCube for the combined data analysis of IC40 and 
IC59. 

The benchmark model assumes constant values for the 
bulk Lorentz factor, the observed variability time and 
the peak luminosity for every burst. Recently, it was 
suggested that there are correlations between the bulk 
Lorentz factor and the isotropic energy, and between 
peak luminosity and the break energy of photon spec- 
trum. Using such inherent relations to derive the Lorentz 
factor and the peak luminosity, we re-calculate the neu- 
trino flux and find that the flux adopting these relations 
exceed the 90% CL upper limit for the assumption of 
t° h = 0.01s for every long burst. This constrains the 
baryon ratio to be r\ v < 10, which, however, could be 
relaxed if the variability times for most GRBs could be 
larger. 

We also calculate the cumulative diffuse flux from 
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GRBs using three different luminosity functions existing 
in the literature. For the luminosity functions of Guetta 
& Piran (2007) and Wandermann'fe Piran (2009), the 
expected flux is below the IceCube upper limit for both 
the case that assumes T = 10 2 5 for every bursts and the 
case considering the inherent relation between T and the 
peak luminosity. However, for the luminosity function 
obtained in Liang et al. (2007), the expected flux exceeds 
the IceCube limit for both cases. The non-detection of 
diffuse neutrinos then constrains the baryon ratio to be 
rj p < 8 in this case. 

GRBs have been proposed to be potential sources for 
UHECRs, besides active galactic nuclei (e.g. Biermann 
& Strittmatter 1987; Takahara 1990; Berezinsky et al. 
2006) and hypernovae/supernovae with relativistic com- 
ponents (Wang et al. 2007; Liu & Wang 2012; Murase 
et al. 2008). Neutrino detection would provide evidence 
for cosmic ray protons in GRBs. On the other hand, the 
non-detection by current IceCube can not yet exclude 
this connectioro, as the required baryon ratio for GRBs 
to be the sources of UHECRs is rj p ~ 5 — 10 (Liu et al. 
2011) for local GRB rate of R ~ lGpc^yi'- 1 (Liang et 
al. 2007; Wanderman & Piran 2009). Future more sensi- 
tive observations by IceCube or other neutrino telescopes 
may put more tight constraints on the baryon ratio and 
would be able to judge the GRB-UHECR connection. 
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Fig. 1. — The neutrino spectrum for one typical GRB obtained, respectively, with the method adopted by the IceCube collaboration 
(Abbasi et al, 2010, 2011a, the IceCube collaboration, 2011) (dark gray solid line), modified Guetta et al. (2004)'s method (blue solid line), 
our modified analytical method (purple solid line) and our numerical method (red solid line). The parameters used in the calculation for 
this GRB are: a = 1, = 2, fluence F° h = 10~ 5 ergcm -2 (in lOkeV to lMeV), z = 2.15, peak energy e° b b = 200keV, peak luminosity 
L 1 = 10 52 ergs — , bulk Lorentz factor Y = 10 2 ' 5 , the observed variability timescale t° h = 0.01s and the baryon ratio r\ v = 10. 



TABLE 1 



L-y (ergs x ) 


r 


z 




10 52 


10 2.5 


2.15 


26.0 






1 


39.9 


L-,G 


r G 


2.15 


8.16 






1 


7.79 


L-/G 


r L 


2.15 


9.07 






1 


7.72 



NOTE— The critical value of the baryon ratio rj p . c for the combined IC40+IC59 analysis obtained by adopting different assumptions for 
the bulk Lorentz factor, peak luminosity and redshift (for long GRBs without measured redshifts). L 7 q represents the peak luminosity 
obtained by using the — L 7 relation in Ghirlanda et al. (2011). Tq and IY are the Lorentz factors obtained with the relations of 
^iso — r in Lv et al. (2011) and Ghirlanda et al. (2011) respectively. Here, the observed variability timescale for long GRBs is assumed to 
be 0.01 s. 
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Fig. 2. — Neutrino spectra calculated numerically by adopting the internal shock radius R = 2F 2 ct° h /(l + z) for 215 GRBs (light 
red lines) observed during the operations in IceCube 40-string and 59-string configurations. We take the same GRB samples, the same 
assumptions for GRB parameters and the same effective area as function of zenith angle as those used by IceCube collaboration. The red 
thick solid line represents the sum of the 215 GRB neutrino spectra and the thick red dashed line is the corresponding 90% CL upper 
limit of IceCube. The thick dark gray solid line and dashed line are the predicted total neutrino spectrum and the corresponding 90 % 
CL upper limit given by the IceCube collaboration for the combined data analysis of IC40 and IC59 (The IceCube collaboration 2011). 
The blue solid and dashed lines correspond to the expected spectra and 90 % CL upper limit obtained by using the modified method in 
Guctta et al. (2004). The purple lines represent our modified analytical calculation as a comparison. For the above calculations, we adopt 
benchmark parameters, such as, the peak luminosity L 1 = 10 52 ergs — 1 , the observed variability timescale t° = 0.01s for long GRBs, the 
Lorentz factor T = 10 2 ' 5 and the baryon ratio r) p = 10 for every GRB. 
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Fig. 3. — The spectra of the total neutrino emission produced by 215 GRBs assuming the same dissipation radius for every GRB at 
_R = 10 12 cm (the black solid line), R = 10 13 cm (the blue solid line), R = 10 14 cm (the green solid line), R = 10 15 cm (the yellow solid 
line),ij = 10 16 cm (the red solid line) respectively. The corresponding upper limits are shown by the dashed lines. Other parameters are 
the same as that used in Figure 2. Note here, the red, green and yellow dashed lines are overlapped with each other, because the spectrum 
shape of the red, green and yellow solid lines are similar in the energy range of 10 5 GeV — 3 X 10 6 GeV. 
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Fig. 4. — The spectra of the total neutrino emission produced by 215 GRBs assuming different fireball parameters in the standard 
internal shock model. The solid red line represents the spectrum that adopts benchmark parameters as in Figure 2. The solid blue line 
represents the spectrum that adopts the relations of Ei so — V and — L 7 in Ghirlanda et al. (2011). The solid green line represents the 
spectrum that adopts Ei so — V relation in Lv et al. (2011) and — L 7 relation in Ghirlanda et al. (2011). The dashed lines are the 
corresponding upper limit by 1C40+IC59. Left panel: z = 2.15 is assumed for long GRBs without measured redshifts; Right panel: 2 = 1 
is assumed for long GRBs without measured redshifts. 
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Fig. 5. — The expected diffuse muon neutrino flux from GRBs and the IceCube limits. The blue, green and red lines represent the 
fluxes obtained with luminosity functions of Liang et al. (2007), Wanderman & Piran (2009) and Guetta & Piran (2007) respectively. 
The solid line and dashed lines correspond to different assumptions about the Lorentz factor used in the calculation, but with the same 
observed variability timescale t° h = 0.01 s for long GRBs and the same baryon loading ratio r] p = 10. The black thick solid line is the IC40 
upper limit on the diffuse muon neutrino flux given in Abbasi et al. (2011b), while the dark gray solid line is the upper limit for one-year 
observation of the complete IceCube, extrapolated from the upper limit of IC40 via 

A IC86 ^ 3Al% 40 (Karle 2011, Hummer et al. 2011). 



